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QED theory of multiphoton cascade transitions in atoms and ions is developed. This theory allows for the 
accurate description of the process important for astrophysical studies of the cosmological hydrogen recombi- 
nation. In particular the 3s — >■ Is + 2y, 4s — )■ Is + 2j and 3p — ^ Is + 87 processes are considered and some 
controversies existing in the literature are resolved. 



I. INTRODUCTION 



The interest to the multiphoton cascade transitions in hydrogen during the last decade was triggered by the accurate mea- 
surements of the asymmetry in the temperature and polarization distribution of the Cosmic Microwave Background (CMB) 1 1 1, 
|2|. The launching of the Planck Surveour enables to perform the measurements with accuracy 0.1%. It is a challenge to the 
theory to perform the calculations of the properties of CMB with the same accuracy. For this purpose the adequate theory of 
the cosmological hydrogen recombination should be developed. The modern theory of this recombination starts from the works 
by Zel'dovich, Kurt and Sunyaev |3| and by Peebles |4|. According to |l3l, IJl the one-photon transitions from the upper levels 
to the lower ones did not permit the hydrogen atom to recombine, i.e. to reach the ground state. Each photon released in the 
one-photon transition in atom was immediately absorbed by another atom. This reabsorption process did not allow the radiations 
to escape the interaction with the matter However if the atom anives in the 2s-state, then it decays via the two-photon transition. 
These two photons escape the reabsorption and the recombination occurs. It was first established in fj|, |4| where 2s — Is 
transition was found to be the main channel for the radiation escape and formation of CMB. Hence the recent properties of the 
CMB are essentially defined by the two-photon processes during the cosmological recombination epoch. 

Apart from 2s — Is transition, as it was noted recently in |5| the two-photon decays from the excited states with the principal 
quantum number n > 2 also can contribute at the 1% level of accuracy. This idea was further developed and intensively 
discussed in ||6|-|[8|. There is a difference between the decay of ns {n > 2), nd states and the decay of 2s state. This difference 
is due to the presence of cascade transitions as the dominant decay channels in case of ns [n > 2), nd levels. For the 2s level the 
cascades are absent. The cascade photons can be effectively reabsorbed and therefore the problem of separation of the "pure" 
two-photon emission from the cascade photons arises in connection with the escape probability. 

A problem of cascade separation appeared to be nontrivial and caused several controversies in the literature. For the first 
time this question was raised in ||9l for the two-photon transitions in the two electron Highly Charged Ions (HCI). The same 
problem was considered later in |10|. In 111], |12| a general QED approach was developed which allowed for the description 
of the few-photon cascade transitions. This approach was based on the F. Low theory of the line profile in QED |13|. The 
new interest to cascade separation problem did arise in the context of the cosmological hydrogen recombination. The ambiguity 
of this separation was shown in [ 14 1 where it was demonstrated that the contributions of the cascade, "pure" two-photon and 
interference terms can vary depending on the method of calculation and only the total transition rate remains invariant. 

The controversies in the cascade description appeared in connection with the cascade regularization methods. It is convenient 
to describe these controversies using the simplest example: the two-photon transitions 3s — > Is + 27. In hydrogen we can 
restrict ourselves with the nonrelativistic theory and only the electric dipole (£'1) transitions. In 3s Is + 27 transition a single 
cascade 3s 2p + 7 — !• Is + 27 should be taken into account. 
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A total transition rate WTj_ig can be written as 







where dWl]_i^{uj) is the differential transition rate, uj is the frequency of one of the emitted photons, = E^g — Eis. The 
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differential transition rate dWl]_ig{ijj) consists of three terms: cascade contribution, "pure" two-photon contribution and the 



interference contribution: 



The cascade contribution can be presented as the sum of the contributions of two cascade links (resonances): 
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where two resonant frequencies are: uj^^'^^ = E^s — £'2p and oj'''^*^ = E2p — Eis- The corresponding resonance contributions 
were presented in lil4J on the basis of the QED approach developed in ifTTI . lfT2l . These contributions look like 

27(resonance 1) _ ^3s + ^2p ^37-2p (^''''''^ )^2p- Is iuj'^'"''^)dui 
„,,2,iresonance 2) _ ^37-2^ K^^^ ) ^^2^-1^ K^^^ )'^^ 

Note that the factor ^^'p^^^'' in Eq. (4l was lost in |14| which led to the wrong values for PFgJ^ig, ^ll-xs transition rates 

different from correct ones in the third digit in case of 3s. This mistake was noticed in ifTSl . 

Here Fs^, are the total widths of the levels 3s, 2p and W\2-ip^ ^2^-is ^^"^ 'he one-photon transition rates. In the 

nonrelativistic limit = W\2-ip-: = W2p-is- Then, integrating Eqs. mI, i5« over uj and taking into account Eq. ill we 
find l-J LJ l_l 
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From Eq. (j9jl follows that the deviation of W^]_i^ from T^s is quite small; actually this deviation arrives only in the 5th digit. 

A controversy did arise in connection with the regularization of the divergent cascade terms in the integral (1). In Q, Q, 
[|8] as in some further papers the regularization was performed by introducing the widths for the intermediate np states, i.e. 
replacing the energy Enp by Enp ~ |r„p. This replacement was made phenomenologically within the Quantum Mechanical 
(QM) description. In case of 3s — 2p — Is cascade T„p — F2p. Therefore the contribution of the resonance 1 instead of Eq. Q 
looked like 

„.r2',{resonance 1) _ ^3s'-2p (^'^'"'^ ) ^2p- 1^ {uj'^'"''^ ) duJ 

while the contribution of dW^]!f2p°i^"'^'^ remained the same as in Eq. p^. In principle. Eq. (|lo|) for the cascade transition 



3s — 2p — Is can be considered as an approximation to Eq. Q since 

r3s + r2p _ , r3s 



r2p(i + T^)^F2p(i + o.oi) . (11) 
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However, S.G. Karshenboim, V.G. Ivanov and J. Chluba flS] insisted that the correct expression for the contribution of the first 
resonance to the W^]_i^ transition rate is Eq. ( 10 1 but not Eq. 

In the present paper we demonstrate explicitly that our expression Q is exact and the expression ( 10 1 can be considered only 
as an approximation. Moreover we prove that the statement made by S.G. Karshenboim et al in Iil5i concerning the independence 
of the two-photon transition rates on the initial state widths, is definitely wrong. 

Note that we do not deny, in principle, the usefulness of QM approach; still this approach should be applied with more caution. 
All these circumstances require full clarification which will be given in the present paper Completing the introduction we have 
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to stress that the insertion of Eq. ( [TO) l in the integral Eq. ([T|l gives exactly the same result Eq. (|6| as the insertion of Eq. 
Q. One could think that both methods of regularizations, the QED one applied in IHI, |[12J, L14|, |16|, |17j and the QM one 
applied in fT], IH, ifTSl . are equivalent. However, this is not the case for three reasons. First, this equivalence for the cascade 
contributions is approximate. Eqs. (j6|, (j?]) are valid up to the small corrections of the order Fss/ojo. With the same accuracy holds 
the mentioned equivalence. The statement that the cascades do not contribute at all to the radiation escape is also approximate 
since the cascades can not be separated exactly from the "pure two-photon" contribution. Second, the interference contribution 
also requires regularization and depends on the regularization method. The integral contributions from the interference terms, 
unlike the integral cascade contributions are not equivalent for the different regularization schemes. Third, in the astrophysical 
applications the frequency distributions for the two-photon decays are converted usually with some other functions. This also 
violates the equivalence mentioned above. The total contribution of all the excited states to the radiation escape according to |6| 
is about 0.4%. Thus the error due to the employment of the wrong regularization scheme hardly can exceed 0.1%. Nevertheless 
having in mind rapidly growing accuracy of the astrophysical measurements of the CMB the development of the accurate QED 
theory of the processes in hydrogen, connected with the cosmological recombination seems to be necessary. 

The paper is organized as follows. In section 11 we start with the QED derivation of the Lorentz profile for the one-photon 
transition from the excited state to the ground state. This derivation repeats shortly the derivations in |lT|, [TZ] but is necessary 
to introduce the basic formula and notations. As an example the Lyman-alpha transition 2p — Is is considered. In section III the 
two-photon transition rate to the ground state from the ns-state in the presence of cascades is described in general. In section IV 
the regularization of the two-photon transition 3s — Is is analysed. For this transition the QED regularization scheme fTTl, fT2l 
deviates from the QM approach employed in |6|-|8|, 1 15|. In section V the "pure" two-photon and interference contribution to 
the 3s — Is two photon transition rate are described. In section VI the same derivations are made for the two-photon 4s — Is 
transition: there is an important difference between 4s — Is and 3s — Is two-photon transitions due to the existence of several 
cascade channels in case of 4s — Is. In section VII the 3-photon transitions are analysed with 3-photon decay 3p — Is as an 
example. Concluding remarks are presented in section VIII. 

II. ONE-PHOTON TRANSITION TO THE GROUND STATE 

The full QED description of any process in an atom should start with the ground state and end up with the ground state too, 
i.e. the excitation of the decaying state should be always included. For the resonant processes, e.g. for the resonant photon 
scattering the absorption part of the process can be well separated from the emission part, so that the description of the decay 
process independent on the excitation becomes possible. In this way the theory of the multiphoton processes in atoms was 
developed in fTTl, fT2l. A simplified version of the theory which starts directly from the excited state was considered in fTSl. 
This version allows for the correct description of the complicated multiphoton processes with cascades but does not allow to 
trace down the details of the regularization of the divergent cascade contributions, i.e. one has to refer to the more elaborate 
evaluations 1111 . |il2|. Since the controversies mentioned above in section I, concern namely the regularization methods, in this 
work we follow the description formulated in ifTTI . lfT2]| . 

Having in mind the recombination processes in hydrogen atom we consider first the resonance photon scattering on the ground 
Is state with resonances corresponding to the np states. In our derivations we will fully neglect the photons other than El which 
is reasonable for the neutral hydrogen. It is important to stress that we consider the free atoms which are excited by the photons 
released by the source which line widths is comparable (or larger) then the natural line widths of the resonance atomic state. Thus 
we exclude the special cases of the excitation by the laser with the narrow bandwidths or something equivalent. Our condition 
(broad source width) should correspond the cosmological recombination situation when every atom is excited by the photons 
emitted by another atom. The Feynman graph corresponding to the resonant photon scattering is depicted in Fig. la. 

The S'-matrix element, corresponding to Fig. la, i.e. second-order scattering process, looks like 

where 

^a{x) ^ Mr^e-'""^' , (13) 

4'Ai'r) is the solution of the Dirac equation for the atomic electron, Ea is the Dirac energy, = ■0^70 is the Dirac conjugated 
wave function, 7^ = (70, 7) are the Dirac matrices and x = {f, t) is the space-time coordinate. The photon field or the photon 
wave function j4p(a;) looks like 

aA?) = y^e^e^(£--— *) ^ y^e-'"*4^-T(r) , (14) 
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where e^t are the components of the photon polarization four- vector (e is 3-dimensional polarization vector for real photons), 
k = {k, u) is the photon momentum four-vector, k is the wave vector, uj — \k\ is the photon frequency. Eq. ( 14 1 corresponds to 



the absorbed photon and A*J'''^' corresponds to the emitted photon. Finally, the electron propagator for the bound electron it is 
convenient to present in the form of the eigenmode decomposition with respect to one-electron eigenstates ifTSl 

1 



^ ' ^ 2TTi J ^ uj-En{l-iO) 

Insertion of the expressions ([T3]l-([T5| into Eq. ( [T2| and integration over time and frequency variables leads to 

r,(2)sc „ .r/ N 2 WM^M Jlsn\l^lA^l )nls 

Sis = -2T:i5{uf - u,)e ^ , , , p ^ ■ 

n ^ 

The amplitude U of the elastic photon scattering is related to the S'-matrix element via ifTSl 



(15) 
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Accordingly, we will obtain the scattering amplitude 



{2)sc 



+ Els — En 



(16) 



(17) 



(18) 



where the energy conservation law implies that |fcy | — 

For the resonant scattering process the photon frequency uji — ujf is close to the energy difference between two atomic levels, 
hi case of np resonance uji ~ E. 



np 



Els- Accordingly we have to retain only one term in the sum over n in Eq. ( 18 1 



'2)sc 
Is(np) 



2 WM^M )lsnp\l^i-^ii 



npls 



LUf + Els — E„ 



(19) 



Eq. ( 19l reveals that in the resonance approximation the scattering amplitude is factorized into an emission and absorption 
parts. The energy denominator should be attached to the emission or absorption part depending on what we want to describe: 
emission or absorption process. In particular, the first-order emission amplitude can be expressed as 



7- rem 
^ npls 



^ \ /M M Jlsnp 
UJf + Els — Enp 



(20) 



The nonresonant corrections to the resonance approximation, first introduced in fT3l were recently investigated in fT9l-f2Tl. 
The role of these corrections appeared to be negligible in most cases. These corrections arise when one takes into account the 



terms other than the resonant one in sum over n in Eq. ( 18 1. The same concerns nonresonant contribution to the scattering 
amplitude which arises when we interchange the position of the photon lines in Fig. 1, i.e. when the emission of the photon 
occurs prior to the absorption. 



The energy conservation law which follows from Eq. ( 17 1 reads 



The resonance condition one can write in the form: 



E, 



np 



Eis\ = \ujf-Enp + Eis\ 



np 



In cases, when we can neglect r„p in Eq. (22 1 this equation takes the form of the energy conservations law 



<^f — Enp — Els 



(21) 



(22) 



(23) 



In particular we can use Eq.( 23 i in the numerator of Eq. ( 20 1 but not in its denominator. 

To derive the Lorentz profile for the emission process we follow the Low procedure fVP\, i.e. insert infinite number of the self- 
energy corrections in the resonance approximation into the electron propagator in Fig. la. The first term of the corresponding 
Feynman graph sequence is depicted in Fig. lb. Employing the photon propagator in the Feynman gauge in the form 



D,. 



,Xxi - X2) 



1 

27h 



-t2) 



(24) 
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I,,,. - ^e^lf^l'-- , (25) 

7-12 



where x = (r, t), ri2 = \ fi — 7^2! and defining the matrix element of the electron self-energy operator as ll22ll 
we obtain the following expression for the correction to the scattering amplitude lfT2l : 

The resonance approximation implies ni ~ n2 — np. Then taking into account Eq. ( [T9] ) we can write 

.r(4)sc. _ rr(2)sc. (S(a;/ + Eis))„p,np ,^„. 
^l.(np) - ^l.(np) + iJ^^ _ £;„p ■ ^^^> 

Repeating these insertions in the resonance approximation leads to a geometric progression. Summation of this progression 
yields 

j-j-{4:)sc. _ 2 7m^A' jlanpl,7M^A' Jrtpls ,^„,^ 



(w/ + Si, - E2p - (S(a;/ + 



The emission amplitude looks like 



+ — — + Eis))np,np 

The operator ^(wy + can be expanded around the value Uf + Eig — E^p 

t{iOf + El,,) = Y:{E„p) + (iOf + Els - Enp)^{Enp) + ... , (31) 

where Yi'(Enp) = ;^S(^)|^=_E„p. The first two terms of the expansion (jsTj) are ultraviolet divergent and require the renormal- 
ization. The methods of the renormolization in the bound electron QED are described, for example in [23 J . In order to obtain 
the line profile for the emission process we retain the first term of the expansion (31 1 and consider the energy denominator in 
Eq. p9]l as a complex quantity: 



(^(-^np))np,np — Ej^p 2"^^P ■ (32) 



Here is the electron self-energy contribution to the electron Lamb shift and r„p is the one-photon radiative level width 
II22I . Apart from L^^ contribution there is also the vacuum polarization L^^ contribution lfT2l . but the vacuum polarization 
contribution is pure real and does not change the imaginary part in Eq. ( 32 1. Now the emission amplitude reads 

K^_ls = e )lsn^^ ^ ^33^ 

LOf+Eig — Enp — Lnp + 2 Trip 

where L„p = L^^ + L^^p, 

As a next step one has to take the amplitude Eq. (33 1 by square modulus, then integrate over the photon emission directions 
i7f and sum over the photon polarizations. Defining the one-photon transition rate for the transition np — Is like 

W^;_is = 2^ul^s J (^l(7.<'^""^' Vpi.P , (34) 



where oj^es is the resonant photon frequency. In Eq. ( 34 1 it is assumed also the summation over the degenerate substates of the 
final state and averaging over the degenerate substates of the initial state. These operations we will not designate explicitly since 
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it does not influence our argumentation. The same wifl concern the two-photon and three-photon transitions in the subsequent 
sections. 



From Eq. (33 1 we obtain for the absolute probability of the photon emission with the frequency in the interval between tjjf 
and ojj + do;/ 

dwnp-is{ujf) —- — — r^^— TF^ • (35) 

2n [ujf + Els - Enp - LnpY + jTi 



Due to the factor ^ the Lorentz profile Eq. ( 35 i is normalized to unity for the Lyman-alpha transition 

dw2p-is = 1 • (36) 



I' 



In case n > 2 the Lorentz profile is normalized to the branching ratio for the transition np — Is: 



dwnp-is = —^^-^ = . (37) 



The Lamb shift for the ground Is state enters the energy denominator in Eq. (35 1 in different way. Insertions of the electron 
self-energy corrections in the outer electron lines in Fig. 1, unlike the insertions in the internal electron line lead to the sin- 
gularities when the intermediate states in propagators are equal to Is. This singularities are not connected with the frequency 
resonances. To regularize these singularities one has to introduce Gel-Mann and Low (24] adiabatic S'-matrix as it was done 
in |25J. It was demonstrated that the summation of the infinite series of the singular in the adiabatic parameter A terms can be 
converted to the exponential factor. The amplitude Eq. ([29]l should be replaced by 



/ /i*(fc/e/)\ / /i(feiei)s 

lim C/LVp)(A) = e2^0VV )isnp[lpA^ ^3^^ 

A-i-O ojj: + Els + Lis — Enp — Lnp + ^^np 



Since for the ground state the matrix element (Y,{Eis))isis is pure real, for the probability this gives 



lim |e"i(^(^i=»i=iH = 1 (39) 

A^O 



and thus the Lamb shift Lis arrives in the expression ( [35| for the Lorentz profile. Note, however, that if we apply Eq. ( |39| to the 
excited state and take into account the width of the excited level we will obtain zero transition probability. Strictly speaking this 
means that it is incorrect to evaluate the transition probabilities via the nondiagonal S'-matrix elements as is usually done in QED 
for atoms and it is necessary to start with the process of excitation using the procedure described in the present paper. However 
in most cases Eq. (|39]) can be ignored and the correct results for transitions rates are obtained in a standard way, evaluating 
the square modules of the nondiagonal S'-matrix elements. Only in the special situations as in case of the multiphoton cascade 
transitions considered in the subsequent sections of the present paper, more refined analysis is required. 

in. TWO-PHOTON ns - Is TRANSITION 

In this section we describe the two-photon transition to the ground state using as an example ns — Is two-photon transitions. 
According to our approach we have to start with the Feynman graph depicted in Fig. 2a. The two-photon resonant excitation is 
the most natural and convenient way to describe the excitation process in this case. The resonance condition is 

Wjl + Wj2 = = Ens - Els ■ (40) 

Constructing the S'-matrix element corresponding to the Feynman graph Fig. 2a, inserting the expressions for the electron and 



photon wave functions as well as the expressions for the electron propagators Eqs. (13 i-( 15 i, integrating over time and frequency 



variables and using Eq. ( 17 i for the scattering amplitude results 



" „i:„A^f2+Eis~En,){uf,+Uf,+Eis~En,) 



(a;/2 + w/i - + Els ~ En^) 
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The energy conservation in this process is implemented by the condition 

^fi + = ^ii + (42) 
and the resonance condition is given by Eq. ( [40| . From Eq. (|40]) follows the approximate energy conservation law similar to Eq. 

\uJf,+UJf^- Ens+Eu\i^rns , (43) 

which can be replaced by equation similar to Eq. ( [23| ) 

uJf,+ujf^^E„,~Eu, (44) 
when r„s can be neglected. According to Eqs. (|40]) and (|42]) the last energy denominator in Eq. can be replaced by 

'^f2+ ^fl - + Els - = UJi^ + Els - Ens i (45) 

i.e. does not depend on the frequencies of emitted photons. 

In the resonance approximation we retain only one term n2 — ns in the sum over n2 which yields 



U. 



(A)sc 
Is(ns) 



(CJ/^ + Els " -B„ J(w/2 + + £^15 - Ens) 



rials 



{uJi^ + Els — En^) 



Starting from Eq. (46 1 we can write down the expression for the two-photon emission amplitude as 



u. 



(2)e 



rii 



(ujf^ + Els - Eni){u}f^ + W/i + £^15 - Ens) 



(46) 



(47) 



with the condition Eq. ( 40 1 remaining valid for the frequencies ujf-^^ujf^ due to Eq. ( 42 1. 

Up to now all formulas above in this section were valid for any ns levels, beginning from n = 2s. Now we have to take into 
account the form of the resonance produced by the second energy denominator in Eq. (47 1. The width of this resonance for 2s 
level is defined by the two-photon transition 2s — > Is + 27. This width should arrive as the imaginary part of the matrix element 
of the second-order electron self-energy operator, i.e. from two-loop insertions to the Feynman graph Fig. 2a. The rigorous QED 
derivation of the two-photon widths from the contributions of the two-loop Feynman graphs is stiU absent but the contribution 
of the 2s level to the CMB history is very well known from [3J, [4| and later works. 

Therefore, we will restrict our studies with n > 2. For n > 2 there is always leading one-photon contribution to the total 
width Tns, for example the W^]_2p transition rate in case n = 3. Assuming the existence of such a contribution we will continue 
our studies by inserting the one-loop electron self-energy corrections to the central propagator in Fig. la (the Low procedure). 
The first term of the Low sequence is depicted in Fig. 2b. 

Returning back to the scattering amplitude Eq. ([45]l and proceeding along the same way as in the case of the one-photon decay 



we obtain an expression similar to Eq. ( 27 1 in the one-photon case: 



(6)sc 
Is(ns) 



E 



lsniW/i2^M2 M1TI2 



{U}2 + Els - )(t^/2 + W/i + Els - En^) 



(48) 



+UJf^ +Eis))r 



{UJI + Els - EnJ 



i^h +^fl +Els " Ens) 

In the resonance approximation setting n2 = = ns we have 

(6)sc _ jj(i)sc (S{a;/i + W/2 + Eis))nsns 



isins) (^^^ +^^^+Ei,-Ens) 

Producing further the Low sequence and performing the summation the arising geometric progression results 



TTS. 



(ns) - e 2^ 

ni 



Els - Ens — {^{'^fi + + Eis))r 



i^f2 



Els — Eni) 



(49) 



(50) 



(Wij + Els — En 
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The emission amplitude for the two-photon decay process ns — ^ Is + 27 looks like 



V f^ti^Ml )nins\ f;j2^M2 )nin2 ^ (51) 



1 



^f2 + ^fl + ^Is — Ens — C^i^fi + ^f2 + Eis))ns.ns 

At the point of the resonance we expand the operator 

^{ujf,+ujf,+Eu)^f^{Ens) + - (52) 

and using the equality 

i'^iEns))ns,ns = L^f — 2^ns , (53) 



arrive at 



jjem g2 V l^^l^^J■l )nins\ Ip,2^l^2 (54) 

i^f2 + Els - En,)iLJf, + ujf, + Els - Ens + ^Tns) ■ 



In Eq. ( 54 1 we have omitted the Lamb shift of the ns level in the second energy denominator In what follows the Lamb shift 
will play no significant role in our derivations. 

The value Tns is defined in a different way for the different ns states. For example, for n = 3 T^g = M^3^7_2p since there are 
no other one-photon decay channels for 3s level. The further investigations of the two-photon transition probabilities should be 
performed separately for different n. In the next section we will continue these investigations for 3s — > Is + 27 transition. 



IV. TWO-PHOTON 3s - Is TRANSITION 



The further studies of the 3s — Is transition we can start with the expression for the emission amplitude Eq. (54i written for 



the case ns — 3s. The Feynman graphs for the resonance two-photon scattering with the excitation of 3s level are depicted 



in Fig. 3. To the expression Eq. (54i we have to add also another term corresponding to the Feynman graph Fig. 3a with the 



interchanged positions of the /c/^ , e/j^ and /c/, , e*/, photons. This yields 



-/2 + Els En, + U,,+Eis-En, ^''^ 



1 



+ "^/i + Els — E^s + IFss 

For the 3s — Is two-photon transition only one cascade is possible: 3s ~ 2p — Is. Accordingly, the two new resonance 
conditions arise (these resonances were defined also in section I): 

^'■"■1 = E^s - E2p , (56) 
c^""-' - E2p - Els ■ (57) 



Consider first cascade contribution to Eq. (55 1. For this purpose we have to set ni = 2p. Then 



jjem., cascade — g2 J Wmi^Pi j ls2p (,7^12 ^M2 )2p3s ^ Ufii^f^i )ls2p[^H2^t^2 )2p3s . ^ ^^g-^ 

3s-2p-ls 1 _|_ j^^^ _ j^^^ j^^^ _ j^^^ 

1 



'^f2 + ^fi + Els - Ess + fFss 
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The first term in tlie curly brackets describes the resonance ( 56 1, the second term describes the resonance ( 57 1 (see Appendix 
A). Applying the Low procedure (insertions and summation of the infinite chain of the electron self-energy corrections) to the 
upper electron propagators in Fig. 3b we find 



u: 



em. casca 
3s-2p-ls 



de 



1^, /l*('=/2^/2)\ /„, 



(k 



2p3s 



(k 



1^, A*i^!2'^S2. 
Isrii V /A'2^M2 



2p3s 



^f2 



E 



Is 



2p 



^fl + ^Is — E2p + |r2p 



(59) 



^/2 + ^/i + Els - E: 



3s 



3s 



Now we take U^^2'p^_^is " square modulus, integrate over the emitted photons directions and sum over the polarizations 
of both photons. The formula (34i should be used for presentation of the results of these integrations and summation via the 
one-photon transition rates. Consider first the square modulus of the first term in the curly brackets and the factor outside the 
curly brackets in Eq. (59i. This term is represented by Fig. 3a and corresponds to the contribution of the resonance 1 in Eq. 
(|56]l. In this case we are interested to derive the Lorentz line profile for the upper link of the cascade 3s — 2p — Is. Therefore 
we have to integrate first over frequency of the second emitted photon, i.e. uif^. In principle the integration over both photon 
frequencies should be done with Eq. (|43]l taken into account, i.e. 



dujf 



1 



dujf 



(60) 



where a;,, 



E, 



Els 



Eq. ( 60 1 holds due to the symmetry of Eq. (59 with respect to permutation w/j ^ ujf^. 

The integration over the frequency uj in Eq. ( 59 1 we perform in the complex plane. Since only the pole terms contribute we 
can extend the interval of integration to (— oo, +00) and not to refer to Eq. (43 1 or (6O1. Then using Cauchy theorem after some 
algebraic transformation (see for details the Appendix A) we obtain the cascade contribution (resonance 1) to the differential 
branching ratio 

1/I7 f, ,res.l\Ti/l7 /, .res.2\ 
^3s-2p[^ l^2p-lt 

(here we have changed the notation for the frequency from w/j to uj). 

The differential branching ratio db"^"^ is connected with the differential transition rate dw^g_ig{uj) via 



db. 



resonance 



27 ( 

3s— 2p— Is 



1) 



1 r, 

2ti 



3 s 



r 



2p 



r3sr2p 



'-)dLO 



(61) 



dw'^ns-ls 

r 

ns 



(62) 



This definition concerns not only the cascade contributions but all the contributions in Eq. ([2]) for the two-photon decay of 
any ns-state: 

1 



db^'^ 



,,2j{cascade) j,2f (pure) 2-y (inter f erence) _ , ^-,y^^seu.u.e, , ^lypui-e, , 

"■"ns-ls "T ""ns-ls "r""ns-ls " -p l""'ns-ls "'"""'ns-ls "'"ns-ls 



27(p«re) 



2j(interf erence) 



(63) 



Combining now the formulas ( [6T| , ( 62 1 we arrive at the expression (j4]l presented in the Introduction. The integration of Eq. 



(63 I over the remaining frequency will give the total branching ratio 



blZ 



27 



Note that this last integration according to Eq. (60 1 should be done within the interval (0, w„ 
mation can be used. 



(64) 

c) since now no pole approxi- 

The second term in the curly brackets in Eq. ( |59] l is represented by the Feynman graph Fig. 3a (with the change of the photons 
u!f-^ ^ ujf^) and corresponds to the resonance 2 in Eq. (57 1, i.e. to the lower link of cascade. To obtain the Lorentz profile for 
this lower link we have to integrate over the frequency of the first emitted photon, i.e. again over ojf^ after taking the square 
modulus of this term and the factor outside the curly brackets in Eq. (59i. Replacing notation ujf-^ to ui we obtain the cascade 
contribution (resonance 2) to the differential branching ratio: 



77 2^(resorLance2) , \ 
dhs~2p-ls (^) 



1 1 



27r r 



3s 



,.es.2)2 ^ lr2^ 



(65) 



Combining the formulas (62 1 and (65 1 we arrive at the expression (jSj presented in the Introduction. 

The interference between two terms in Eq. (59 1 should not be taken into account since these two terms correspond to the 
resonances located far from each other: at the distance ^max in the frequency scale. 
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V. "PURE TWO-PHOTON" AND INTERFERENCE CONTRIBUTIONS 



Returning to Eq. ( 55 i we consider this expression with the state 2p excluded from the summation over ni as a "pure two- 
photon" contribution to the transition amplitude 3s — Is. This corresponds to the "pole approximation" employed in Section IV 
for the description of the cascade contribution: extension of the of the first frequency integration over the interval (— oo, oo). In 
llT4l the more general approach was developed, when the resonances were regularized only within the "windows" of the different 
breadth. Then the 2p state should be eliminated from the sum over ni only within "windows". The "pole approximation" 
corresponds to the window breadth [cj] = oo. This case was considered in [SJ. 



Since the energy denominators (apart from the factor outside the curly brackets in ( 55 1) now become nonsingular we can 



employ the energy conservation law Eq. (42 1 to replace the frequency cj/^ in the second denominator in curly brackets in Eq. 
(55 I by cj/j = i^max ~ ^f2 - Then the "pure two-photon" contribution to the amplitude becomes 



em.^pure 2 



E 



17^1^^! jlsnAl H2^tJ.2 jni3s ITah^Mi j Isrii 17ai2 ^^2 jni3s 



^f2 + ^Is ~ 



1 



X (66) 



^f2 + ^fi + Els — E^s + fTss 



Taking Eq. ( 66 1 by square modulus, integrating over the directions of the emitted photons, summing over the polarizations, 
integrating over uj , and changing the notation ~ lo results 



E 



ITmi^Mi nsniW/J2^M2 )ni3s [lfJ.i^fJ.1 JlsriiWM2^M2 Ml3s 



e/l e/2 
,*(fe/2e/2). 



duf^ dvf^ 
(27r)3 (27r)3 



+ Els — -E'm 



E^s-i^- E, 



ni 



-duj 



3s 



(67) 



Then, according to Eq. (63 i the "pure two-photon" contribution to the differential transition rate is 



dW!]lTs-\^) - -^)^EE / X 



E 



Vly-l-^V^l jlsni (,7/J2^M2 )nx3s ^ (.7^1 ^M l j Isni (,7^12 ^/'2 jni3s 



W + Els — En^ 



Ess — oj — E„i 



duj 



(68) 



Now, using Eq. (59 1 and Eq. (66 1 we can write down the interference contribution to the differential branching ratio as 

27(mter/ei-ence) _ r, n [ dvj-^ f dvj^ f 2 , ,2 TTem{pure)* j-j-em,.{cascade) _ 

dbss-is -^^^2^2^Jj^J]2nfJ'^'^fi'^fi'^foy3s-is U ' 



(69) 



E 



Upi^Ml jlsm (7/i2^M2 jrii3s Wah^Mi ) Isrii (7ai2 ^^2 jni3s 



^f2 + Els — En I 



Ess — Eni — W/2 



Umi^A'i jlsnii7M2^M2 Jni3a ^ (.7mi^Mi j Ism (,7At2 ^M2 )»i3s ^ 



^f2 + Els — E2p + |r2p 



+ -Els — E2p + |r2p 
da;/2 



The integration in the complex plane can be extended over the entire interval — go ^LOf^ ^ +00 since only the pole term 
contributes; then we have to take the real part of the expression obtained. 
The "pure two-photon" amplitude in Eq. ([69]l we can assume to be pure real. 
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Then using Eq. (62 1 for the interference contribution to the differential two-photon transition rate 3s — Is we find (changing 
notation oj to a;) 

aH/g^-is - ^ duj + div . (/U) 



j3s _ 2,' 3s \2 J WMl^Ml )lsni\JiJ,2^ti2 Ml3s W/Ji^Ml ; V /M2 ^A»2 ;ni35- 1 ,_^s, 



Vl — 17mi^Mi jls2j9i7M2^M2 j2p3s , (/j) 



V-l — 17mi^Mi jls2pl7Ai2^M2 j2p3s ■ yl"') 



In Eqs. (72 1 - (74 1 it is assumed that \kf^\ = w, |fc/J = [uSq' — w). Unlike the cascade contribution, the dependence on 
+ in Eq. ( 70 1 is essential and does not disappear after the insertion of Eq. ( 70 1 in Eq. (|T]i. The employment of the correct 



expression Eq. (70i for the interference contribution is important for the evaluating of the CMB properties in astrophysics. 



VI. TWO-PHOTON 4s - Is TRANSITION 

Repeating the derivations for the 3s — Is transition for the case of 4s — Is transition we present first the contributions for the 
two cascades: 4s — 2p — Is and 4s — 3p — Is: 

1) Contribution from the upper link 4s — 2p of the cascade 4s — 2p — Is: 

27(resonancel) _ 1 + V 2'p W^47-2p (^''''^'"^ )^2p-ls (^'''"''^)^'^ 



where 



2) Contribution of the lower link 2p — Is of the cascade 4s — 2p — Is 



E^s — E2p, w'''^*'^ — E2p — Els ■ (76) 



3) Contribution from the upper link 4s — 3p of the cascade 4s — 3p — Is: 

dH/4._3p-ls (^_^re..3)2+l(r,^+r3p)2 ' ^^'^ 

where 

_ -g^^ _ JT^^^^ o;'"'^* '* = E^p — Eig ■ (79) 

4) Contribution of the lower link 3p — Is of the cascade 4s — 3p — Is 



.UJ 



dl^4.-3p-l. = _ ^.e..4)2 + lr2^ " ^^^^ 
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Insertion of Eqs. (|75]l, (|77]i, (|78]l, (|80|l in Eq. ([T]i yields 



2^{cascade) 



4s-ls 



"0 4 



''-'' wl:_,^ = + bl;_,,wl]_,^ , (81) 



r 



3p 



where wq = E^s ~ Eis, b^p^is is the branching ratio for the transition 3p — Is. We took into account that T^p — Wg^.^^ 
W':_,^ and blZ. 



3p-2s 3p-ls W^^_,^ + Wi;_,^ 

Hence 



where 



2--) {cascade) 



4s — Is 



4s 



r4. = wl]_^p + M/ilgp 



(82) 



(83) 



rate looks similar to the 3s — Is case (see Eq. (68 i): 



unlike Eq. ([8| in case of 3s — Is transition. The "pure two-photon" contribution to the 4s — Is two-photon differential decay 

(84) 



(27r)3 (27r)3 



E 



V lpi^fJ.1 Jlsni \ lfj.2 



AsniWM2^M2 Mi4s W/Ji^Ml Asni 1/^2 /^2 Mi4s 



E4s — En^ — UJ 



duj 



Interference contribution for the 4s — Is transition consists of 4 terms, corresponding to the four resonances, presented by Eqs. 
(|76ll, (|79l) 



where 



dW, 



2'y {inter f erence) 



-duj 



2{lo -uj''''^^)F^%uj) ^ ^ 2(w-w'''=''4)i^4s(-^) 



e/i e/2 



(85) 



(86) 



ni^2p,3p 



WA'i^Ml Asrai 1.7^2 ^M2 )ni4s 

UJ + Els — Eni 

(Tmi^Mi '^^ '^^ )lsrti (7jJ2^A'2 '^^ )niis 
Eis — En^ — UJ 



(87) 



(pt^ = (7Mi<^'^'''^^)i.2p(7.2^;r'^^'^02p4 



(88) 



^2^ - (7Mi<''^'''^^)l.2p(7.242''^''^^)2p4. 



(89) 



= (7..A;^'^'^^)i.3p(7.242'^^'^^^)3p4. , 



(90) 



= (7m,<''^'''^^)u3p(7m2<2''^''^^)3p4. 



(91) 
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VII. THREE-PHOTON 3p ^ Is + 87 TRANSITION 



In this section we consider the 3-photon transitions, taking as an example 3p ^ Is + 87 transition. The main decay channels 
of the 3p level are 3p — > Is + 7 and 3p — > 2s + 7. Therefore as a resonance scattering process in this case we can choose the 
one-photon absorption and three-photon emission process depicted in Fig. 4. For the transition 3p — >^ Is + 87 there are two 
cascades, containing two-photon links: 8p — > 2p + 27 ^ Is + 87 and 8p 2s + 7 — )■ Is + 87. 

In case of the 8-photon transition we will take into account only cascade contribution. This cascade contribution contains 
necessarily one "pure two-photon" link (3p — 2p or 2s — Is) and is, therefore of the same order of magnitude as the "pure 
two-photon" contribution to the 8s — Is transition. The "pure 3-photon" transitions and the corresponding interference terms 
are essentially smaller than cascade contributions and can be neglected in the "two-photon" approximation |16|. In this sense 
the situation differs from the situation in 3s — Is two-photon decay when we were interested in the "pure two-photon" and 
interference terms. 

The derivation similar to the two-photon case gives the following expression for the 3-photon emission amplitude 3p — Is in 
the resonance approximation: 

(ii;i.-ii;„,+a./3)(i5is-K2+^/3+^/2) 

1 



Els — + LOf,^ + L0f2 + LOf-^ + fFap 



In Eq. ( 92 1 we summed already all the self-energy insertions in the lower electron propagator as is shown in Fig. 4. An exact 
energy conservation law in case of the process Fig. 4 is 

= ujf, + ujf^ + . (93) 

The resonance condition and the approximate energy conservation law in case of 3-photon decay looks like 

\uj, - {Esp - Eu)\ = \ujf, +ujf, + LOf, - {Esp ~ Eu)\ T^p . (94) 
To fix the cascade 3p — >^ 2s + 7 Is + 87 contribution we set n2 — 2s in Eq. ( [92| . This results 

jjem.S-r _ g3 ^ (7^1 ^Ml"^^ ^^)3p2s ^ ^g^^ 

^^"^'^ Els - E-ip + w/3 + LOf^ + ujf^ + IFsp (Els - E2s + ujf.^ + ujf^ + ^T2s) 



(7/J2^A'2^^ )2sni (7/J2^M3^^ 



^ (F F ) +{perm) 



In Eq. (95 I we should include the contribution of the Feynman graphs with all the permutations of the photon lines (perm). 

Now we have to take the right-hand side of Eq. ( 95 1 by square modulus, to integrate over the emitted photon directions and to 
sum over the photon polarizations. Then we have to integrate over the photon frequencies Wf-^, ujf^, Wf^ taking into account the 
condition Eq. (94 1. However when we integrate the contribution of the cascade 8p — 2s — Is we have to take into account that 
the frequency u! is fixed by the resonance condition 

1^/1 - (^3p-^2s)K rsp + Fa. . (96) 

Inserting Eq. ( [96] l in the approximate conservation law (|94]| we obtain the approximate equality 

uJf,+ujf,^E2s-Eis. (97) 

The integration over uj^, ojj^ should be performed in the following way 

d^h j dujf, ^ J dujf., J dcjf^ . (98) 
00 



Eq. ( 98 1 holds after symmetrization of Eq. ( 96 1 via the permutation of the photons with the frequencies ujf^,ujf^. The integration 
over two frequencies (e.g. over ujj^ and in Eq. (95l) can be always extended to the interval [—00, +00] since the pole 
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approximation can be used. The third integration over lo f.^ in Eq. ( 95 1 according to Eq. ( 98 i should be performed over the finite 
interval [0, E2s ~ Eis\. The integration yields 



l37 
"3p- 



.i,(3p-2s-ls) 



dvf 



|jjem.l7(3p_2s-ls) 



(27r)3 (iJgp - E; 



2s 



i(r: 



X (99) 



(27r)3 (27r)3 



E 



WP2^M2 Asrai l7At3^/^3 )ni2s 

Els - E„i + UJf^ 



3p 



2s) 



27 



2s-ls 



r,„r 



3p-L 2s 



where Lo^ax = £'2s - E^. 

The physical sense of the dimensionless quantity 63p_is(3p — 2s — Is) should be discussed specially. This quantity should 
define the 87 transition rate 3p — Is via the channel 3p — ?► 2s + 7 Is + 87. This transition rate is very small compared to the 
main decay channel for the 3p state, i.e. W^p_ig. We assume that the quantity fe3p_is(3p — 2s — Is) is the branching ratio for 

the 3-photon transition rate T4^3p_i5(3p — 2s — Is) to the direct two-photon transition rate W'^'^ = T2s- Then from Eq. (99 1 it 



follows 



<1 



i,(3p-2s-ls) 



3p-2« w27 

*^'^2s-ls 

i 3p 



(100) 



In the same way the contribution of the 3-photon cascade 3p — 2p + 27 — > Is + 87 can be analysed. The final result looks like 

(101) 



637(3p _ 2p - Is) = .^.^I^W?^ 



r2pr3p 



3p-2p ■ 



Unlike b^''{3p — 2s — Is) the quantity Eq. 



101 



should be considered as the branching ratio of the transition rate via channel 



3p — 2p — Is to the total width of the 3p level, i.e. T^p. Then 

W|;_i,(3p-2p-ls) 
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'''^3p-2p 

2p 



(102) 



In the equation in |fT6l . corresponding to the Eq. ( 102 1, the factor W^2pLis/r2p — 1 was omitted. Here we keep it to demonstrate 



that the transition channel 3p 
3-photon) is 



2p — Is is a 3-pnoton channel. Total expression for the transition rate 3p — Is (1-photon and 



M/(37) ^w(i-);) 

* ^ 3p—ls * 3p,ls 



w. 



(17) 



2p,ls rr^(27) 

"3p,2p "I" 

^ 2p 



w. 



(I7) 



3p,2s ^(27) 



(103) 



3p 



This expression coincides with one derived in (16) up to the coefficients before the second and third terms in the right-hand side 
of Eq. ( 104 1. In [ 16| this coefficients were evaluated incorrectly and were equal 3/4. This error was notices by the authors of 
ifTSl . However an expression for M^sp-is given in 1 15 1 is different form Eq. ( 103| l. The equation in ifTSl reads 



i^3p-is = r^^ 



3p— Is 



-27 
3p-2p 



pl7 

^ 3p-2s > 



(104) 



where VlZ^,^, T^^.^, 



3p- 

W^J_r,„, rlZ--}^ — W^J_n^, Eq. (|104| differs form Eq. (|103|l by the last term in the right hand side. Eq. (|104|) represents 



and r^2_o„ are the partial widths for different decay channels. Since T^"' 



3p-2p 



3p— Is 



3p-2p 



3p— 2p' 3p— 2s 



3p-2s' 



the total width of the 3p level. As in cases 3s — Is, 4s — Is decays W^3p_is differs from the r3p (W3p_is 7^ ^sp)- In case 
of 3s level relative difference between W3s_is and r3s was quite small: about 10^*'. This difference became of order of 1 for 
4s — Is. For 3p — Is decay it is again of the order of 1, though the values of the last term in Eqs. ( 103 1 and ( 104 1 differ by many 
orders of magnitude. Within our treatment Eq. ( 104 1 arises when we consider Eq. ( 99 1 as a branching ratio of the decay channel 
3p — > 2s + 7 to the total width T^p. Then, multiplying Eq. (99 1 by T^p and setting Vl^27-is/1^2s = 1 we arrive at 



3p 



-pl7 I p27 

^ 3p-ls I ^ 3p-2s 



(105) 
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Vni. CONCLUSIONS 

In this paper we analyzed the problem of the multiphoton transitions with cascades taking as an example the two-photon 
3s — > Is + 27, 4s — > Is + 27 transitions and the three-photon 3p ^ Is + 87 transition. We proved that the regularization of 
the singularities in the expressions for the cascade contributions to the transition rates includes the widths of the both initial and 
intermediate states. This may be important for the the astrophysical studies of the cosmological recombination and, consequently 
to the connection between the process of the radiation escape from the matter in the recombination epoch and the recent studies 
of the properties of CMB. 



It should be stressed that according to our Eq. ( 39 1, in principle, the absolute transition probability from any excited state to 
the ground state cannot be evaluated with the help of the nondiagonal S'-matrix elements, since the transition probability turns 
to zero. However it does not mean that the standard procedure of the evaluation of transition rates via the nondiagonal ^-matrix 
elements is inapplicable. The standard procedure works in aU cases when it is not necessary to take into account the total width 
of the decaying state. When the total width of the decaying state becomes important we need to employ the Low procedure as 
described in the present paper. 

The necessity of this procedure arises when we need to describe the Lorentz profiles for the decay processes of the excited 
states or when we want to regularize properly the expressions for the multiphoton transition probabilities with cascades. In 
principle, it should be possible to derive the correct results within the phenomenological QM approach as well. However, for 
this purpose one should follow the procedure described, for example, in |27| for the one-photon decay of the level, possessing 
the width: the solution of the non-stationary Schrodinger equation with appropriate initial conditions. This derivation should be 
extended to the case of the two-photon decay. 
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Appendix A: Derivation of Eqs. j61) , ( (65) 



Taking the first term in the curly brackets in Eq. (59 1 together with the factor outside the brackets by square modulus. 



integrating over the emitted photons directions and introducing the shorthand notations 

E2p - Eu = w'""' = AEa , (A.l) 

- Eu = ujo = AEb , (A.2) 

we define the double differential branching ratio as 



Using Cauchy theorem we integrate over uif^ in the lower half-plane where the poles are: 

J/J = AEa - lf2p , (A.4) 

Jf^ = AEb-uj^^-I^^,. (A.5) 

The integration results: 
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db 



2-y{resonance 1) 
'3s-2p-ls 



1 



(A.6) 



V2p[i^Eb - aea -ujf, + |(r2p - T3s)][aeb - aea + |(r2p + rg.)] 



+ 



T3s[AEa - AEb +LOf^- ^{T2p - T3,)][AEa - AEb +^{, + |(r2p + Tg,)] 



duji 



Algebraic transformations then lead to 



T3s\AE 



B 



AE/ 



3s-2p-ls 

i(r2p + T-is)] + T2p[aeb - aea -ujf, + f (r2p + r^s) 



[AEb - AEa - 



Kr2p - TsM^Eb - AEa -iOf,r + lir^p + r3,)2] 

_ 1 w/3;-2pK^'^-')w^2;-i.K^^-') 



27r 



r2pr3s 



(Tss + r2p)[AEB - AEa - w/i + i(r2p - Ta.)] 



[A^B - AEa 



Kr2p - T3s)][{aeb - AEa - + i(r2p + r3.)2] 



(A.7) 



After the cancellation of the factor [AEb — AEa ^ ^fi + | (r2p — ^3s)] in the numerator and the denominator of Eq. (A.7) we 
arrive of the expression (61 1 given in the text. 



To obtain Eq. ( 65 i we have to take the second term in the curly brackets in Eq. ( 59 1, to integrate it over the directions of the 
emitted photons and to sum over the emitted photons polarizations. This gives 



db. 



2-y [resonance 2) 
'3s—2p—ls 



27T [{AEa - L.f, Y + iT%][{AEB ~ - u^f^Y + |r 



i-r2 1 

3sJ 



(A.8) 



We again integrate over loj^ in the complex plane obtaining immediately the result Eq. ([65]l. 
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FIG. 1: Feynman graph describing the resonant photon scattering on the ground state of hydrogen atom. In FIG. la the basic process of the 
resonant scattering with the excitation of np state is depicted. In FIG. lb the electron self-energy insertion in the propagator is made. The 
double solid lines denote the electron in the field of the nucleus (Furry picture of QED), the wavy lines denote the absorbed, emitted and virtual 
photons. 
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FIG. 2: Feynman graph describing the two-photon resonant scattering on the ground state of hydrogen atom, with the excitation of ns (n > 2) 
state and the resonance condition a;i + cj2 = E„b — Eis- In FIG. 2a the basic process of the resonant scattering with the excitation of the ns 
state is depicted. In FIG. 2b the electron self-energy insertion in the central electron propagator is made. The notations are the same as in FIG. 

1 
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FIG. 3: Feynman graph describing the two-photon resonance scattering on the ground state of hydrogen atom, with the excitation of 3s' state 
(resonance condition uji^ + LUi^ = E;is — Ejs and the decay cascade resonances uj^'^"'^ = — E'zp, uf""'^ = E'z-p — -Bis- In FIG. 3a 
the basic process of the resonant scattering with the excitation of the 3s level and decay 3s — 2p — Is is depicted. In FIG. 36 the electron 
self-energy insertion in the upper electron propagator is made. Notation 3~s means that the Low procedure is already performed for this electron 
Une. The other notations are the same as FIGS. 1, 2 
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FIG. 4: Feynman graph describing the resonant process with one-photon absorption and three-photon emission for the ground state of the 
hydrogen atom. In FIG. 4a the basic process is shown and in FIG. 4b the electron self-energy insertion in the central electron propagator is 
made. The notations are the same as FIG. 1 



